Effects of Disorder on Coexistence and Competition between Superconducting and 

Insulating States 
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We study effects of nonmagnetic impurities on the competition between the superconducting and 
electron-hole pairing. We show that disorder can result in coexistence of these two types of ordering 
in a uniform state, even when in clean materials they are mutually exclusive. 

PACS numbers: 74.25.Dw, 71.45.Lr, 74.62.Dh 



I. INTRODUCTION 

At low temperatures, many metals undergo a transi- 
tion into a state with a gap in the single-electron exci- 
tation spectrum and become either superconductors or 
insulators with a periodic modulation of the electron 
charge or spin density. The insulating and supercon- 
ducting (SC) orders inhibit each other by reducing the 
fraction of the Fermi surface available for the gap of the 
competing phase. The balance between the two phases 
is very sensitive to the Fermi surface shape and can be 
changed by, e.g., pressure, doping or magnetic field. 

Nonetheless, in a surprisingly large number of mate- 
rials the SC and insulating states coexist^ This pa- 
per is focused on the coexistence of superconducting 
and charge density wave (CDW) states, observed in, 
e.g., layered transition metal dichalcogenides 2iJ-NbSe 2 
and 2£/-TaS2r i ^ the quasi-one-dimensional compound 
NbSe3j& tungsten bronzes A^WO^t^^ and quarter-filled 
organic materials One of the best studied and best 
characterized CDW superconductors is the transition 
metal dichalcogenide 2#-NbSe 2 . At T d = 33.5 K 
this compound undergoes a second-order phase transi- 
tion to an incommensurate CDW state^i*^ which is 
likely driven by the nesting of a part of the Fermi 
surfaceiiMi The resistivity, however, remains metallic- 
like down to T c = 7.2 K, at which this material becomes 
superconducting^^ and the superconductivity coexists 
with the CDW modulations^ The coupling between the 
CDW and SC order parameters, resulting from the com- 
petition between these two states, was observed in a num- 
ber of experiments. Thus, the suppression of the charge 
density modulation by pressure and hydrogen interca- 
lation results in an increase of T c AiSiii A similar in- 
terplay between the CDW and SC states upon applied 
pressure and doping is observed in NbSe3 and tungsten 
bronzesiSi^ 

In this paper we adopt a rather general, though sim- 
plified, viewpoint on the interplay between the SC and 
CDW states. We assume that it originates from the com- 
petition between two different Fermi surface instabilities: 
the instability towards the electron pairing, which gives 



rise to superconductivity, and the instability towards the 
electron-hole (or excitonic) pairing. Here, we focus pri- 
marily on the effects of quenched disorder on this com- 
petition. We show that even in 'the worst case scenario', 
when the two states compete over the whole Fermi sur- 
face and therefore, in absence of disorder, are mutually 
exclusive, disorder stabilizes a uniform state, in which 
superconducting and insulating order parameters coex- 
ist. While having no effect on the superconducting phase, 
nonmagnetic disorder tends to close the CDW gap before 
completely suppressing the corresponding order parame- 
ter. Disorder induces low-energy states by breaking some 
of the electron-hole pairs. The released electrons and 
holes can subsequently form Cooper pairs, resulting in 
the coexistence of the two phases. 

While in usual s-wave superconductors, non- 
magnetic impurities have little effect on the transition 
temperature^ experiments on electron irradiated 
transition metal dichalcogenides have shown strong 
dependence of T c on the concentration of defects^ 
This was attributed to the interplay between the SC 
and CDW orderings: Similarly to effect of pressure?^ 
disorder strongly suppresses the CDW state, which 
results in the observed increase of the SC critical 
temperature. Theoretically, the combined effect of the 
CDW modulation and disorder on the pairing instability 
have been studied in Ref. j2jj, where an increase of T c 
was found. However, in that paper the amplitude of 
the CDW modulation was assumed to be fixed, which 
is clearly insufficient in view of the strong suppression 
of the CDW state by disorder. In this paper we solve 
self-consistency equations for both SC and CDW order 
parameters, which allows us to study the interplay 
between these two different orders and obtain the 
temperature versus disorder phase diagram of CDW 
superconductors. 

The remainder of the paper is organized as follows: In 
Sec. [n] we formulate an effective model describing the in- 
terplay between the superconducting and excitonic pair- 
ing. The self-consistency equations for the two order pa- 
rameters are derived in Sec. IIHI and in Sec. IIVI we ana- 
lyze the phase diagram of the model. In Sec0we discuss 
the electron-hole symmetry underlying the model and its 
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consequences for the phase diagram. Finally, we conclude 
in Sec. IV1I The details of the derivation of the effective 
model can be found in Appendix A. 

II. THE MODEL 

In the following we consider the microscopic Hamilto- 
nian 

^ = | £ £ 4r fe(-Wx) + C/(X)] 1pj(j 

+S2^^LVW^L'VW > (1) 

CTfj' ) 

describing two types of fermions, one with hole- like dis- 
persion (a-electrons) and another with electron-like dis- 
persion (6-electrons) , £j= ai b(k) = ±(k F — k 2 )/2m (h = 
1), where /i = k F /2m denotes the chemical potential (see 
Fig. Here, in comparison with the models generally 
used to represent CDW systems, two nested parts of a 
single Fermi surface are replaced by two spherical Fermi 
surfaces matching at the Fermi wave vector fci?. The 
excitonic insulator (EI) state is the condensate of pairs 
formed by 6-electrons and a-holes (or vice versa) with the 
zero total momentum^ It is an analogue of the conden- 
sate of electron- hole pairs with the total momentum fi,Q, 
where Q is a nesting wave vector, appearing in the CDW 
state. 

The disorder potential C/(x) encapsulates the effect of 
non-magnetic impurities in the system. Here, we assume 
that the latter is drawn at random from a Gaussian dis- 
tribution with zero mean and variance given by 

(P(x)[f(y))=-^i(x-y), (2) 

where T is the inverse scattering time and vp = mkF /2ir 2 
is the density of states at the Fermi energy. For simplicity 
we have assumed the electron and hole effective masses 
to be equal. 

The interaction term characterized by the coupling 
strength g\ describes the attraction between electrons of 
the same type (e.g. due to the phonon-exchange), while 
the (72 term describes the (Coulomb) repulsion between 
the a- and ^-electrons (31,32 > 0). The attraction be- 
tween electrons favors s-wave superconductivity, while 
the second interaction leads to an attraction between 
electrons and holes and vice versa, favoring the EI state. 
Here, we neglect the inter-band electron transitions due 
to scattering off impurities and electron-electron interac- 
tions, so that the numbers of the a- and ^-electrons are 
separately conserved and fixed by the chemical poten- 
tial. Such terms will formally destroy long-range order 
of the EI phase, corresponding to the suppression of the 




FIG. 1: Schematic picture of the hole-like (a) and electron-like 
(b) bands. 



long-ranged CDW order, due to the pinning of the CDW 
phase by randomly distributed impurities. However, for 
the essentially short-length-scale physics we shall discuss, 
these effects may be neglected. In the absence of disor- 
der, the same model Q has been employed to study the 
competition between SC and EI states for an arbitrary 
ratio of electron and hole densities^ The effect of dis- 
order on the EI state alone has been considered in the 
seminal work of Ref . |2^| , where the analogy of the prob- 
lem to an s-wave superconductor in presence of magnetic 
impurities 2 - was drawn. 

The particular fermion-fermion interactions considered 
in Q - attraction between electrons of the same type and 
repulsion between the a- and 6-electrons - open the possi- 
bility to have simultaneously both superconducting and 
insulating instabilities. A more realistic starting point 
would be a model with attractive phonon-mediated in- 
teractions and Coulomb repulsion between all types of 
electrons. However, it is possible to demonstrate that, 
since the former are retarded, while the latter is prac- 
tically instantaneous, the SC and EI order parameters 
turn out to have a very different dependence on the Mat- 
subara frequency uj. This is clearly shown in Fig. [21 The 
SC order parameter is large at small frequencies, while 
at higher values, it decreases in magnitude and finally 
changes sign when lu is of the order of the phonon fre- 
quency FIq~ By contrast, the EI parameter is large at 
high frequencies and has a dip for |a->| < f2o- I n other 
words, the difference in frequency scales of the attrac- 
tive and repulsive interactions allows both instabilities 
to be present simultaneously. Furthermore, in the weak 
coupling limit and for a weak disorder, i.e. r < fioj the 
frequency dependence of the two order parameters can be 
found separately for u> ~ T and lo > D,q. Furthermore, 
it can be shown that the self-consistency equations for 
the order parameters at to — coincide with the ones 
obtained from the model which therefore can be in- 
terpreted as an effective interacting model. Technical de- 
tails together with the frequency dependence of the two 
order parameters and the explicit expressions for the cou- 
pling constants g\ and gi in terms of the Coulomb and 
electron-phonon couplings can be found in appendix lAl 
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FIG. 2: The dependence of the superconducting (solid line) 
and excitonic (dashed line) order parameters on the Matsub- 
ara frequency ui (see appendix^for details). 



III. ORDER PARAMETERS AND 
SELF-CONSISTENCY EQUATIONS 

Four order parameters describing the SC and EI states 
can be introduced by means of the following anomalous 
averages 



Ala = 9l{tpa]tpal) 

A 2T = -gziipLAt) 



Ait = giiipbiM 



V2j 



Since the numbers of the a- and 6-fermions are separately 
conserved, for homogeneous states, we use the global 
gauge transformation 



to make the SC order parameters, Ai Q and A l0 , real and 
positive. Moreover, as electrons and holes are character- 
ized by the same dispersion, we can require, without a 
loss of generality, that 

A la = Ai 6 = Ai > . 

In case of a spin-independent interaction, as in singlet 
and triplet exciton pairs are degenerate in energy. This 
gives rise to a large symmetry class of transformations 
for the EI order parameter A2 CT . In reality, however, this 
degeneracy is lifted by Coulomb exchange interactions 
and inter-band transitions. Therefore, we will assume the 
exciton pairs have zero total spin, i.e. A21 = A21 = A2. 

Finally we note that, when Ai, A2 ^ and A2 has an 
imaginary part, a pairing of electrons of different types, 
Aiab = ~g2(^a-\i>bi) = -g2(ipb-\ipal), may be present^ 
However, one can show the energy of the state with co- 
existing SC and EI orders to be the lowest for real A 2 , 
in which case Ai ao = 0. 



By analogy with the case of magnetic impurities in s- 
wave superconductors^ restricting attention to the limit 
in which the disorder potential imposes only a weak per- 
turbation on the electronic degrees of freedom (/1 <C T) , 
the mean field (saddle-point) equations together with the 
self-consistency equations for the EI and SC order pa- 
rameters can be obtained using the diagrammatic tech- 
nique. However, we will find it more convenient to use 
a path-integral approach. This will also allows us to ob- 
tain straightforwardly an expression for the average free 
energy. 

The quantum partition function, Z — tr[e~^ w ], where 
(3 = 1/T, can be expressed as a coherent state path inte- 
gral over fermionic fields. In order to facilitate the aver- 
aging of the free energy over the disorder potential J2J, 
it is convenient engage the replica tricki 2 ^ 



F = -i(l„2) = 



P n->0 



Once replicated, a Hubbard-Stratonovich transformation 
can be applied to decouple the interaction terms in the 
Hamiltonian. As a result, one obtains: 



Z n = J Z)(*,* t ) J DA 1 DA 2 

{ rP r r , , a 2 a 2 

exp I J dr J dx *t [d T + H J W + 2-i + 2-^ 



Here, omitting the replica indices for clarity, the 
fermion field is arranged in a Nambu-like spinor, ^ T — 
(tpbt, ^a^^l^^li), in such a way the single quasi- 
particlc Hamiltonian takes the following form: 

H = ipT 3 a 3 + C/(x)t 3 + Am + A 2 T 3 cri , (3) 

where, £p = -Vj/2m — \x and the Pauli matrices r c and 
a c (c = 1,2, 3) act, respectively, in the particle-hole and 
the 6, a subspace. 

The ensemble average over the quenched random po- 
tential distribution J5J induces a time non-local quartic 

interactions, J dx(J Q dr^r^) 2 , which can be decou- 
pled by means of a Hubbard-Stratonovich transforma- 
tion with the introduction of a matrix field, S(x), local 
in real space, and carrying replica, Matsubara (uo n = 
(2n+l)ir/f3) and internal (particle-hole and 6, a) indices. 
Integrating over the Fermionic fields 'J, one obtains the 
ensemble averaged replicated partition function: 



(Z n ) 



DA 1 DA 2 / DYjc 



-f3F 



where F is the free energy of the system 

[P f ( A? Al 
0F = dr dx[ 2—± + 2—±- 
Jo J \ 9i 92 



tr In 



(-O- 1 )-^- /dxtr[S(x)r 3 ] 2 (4) 
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and Q is the quasi-particle matrix Green function in the 
presence of disorder: 

-G~ x = -iu n + S,pT 3 a 3 + Am + A 2 T 3 cri + E(x) . (5) 

The matrix field E(x) represents the contribution of the 
non-magnetic impurity interaction to the self-energy. 

The saddle-point associated with the action Q ob- 
tained by variation with respect to the self-energy E, 



E(x) 



1~KVF 



T 3 (x|£|x)t 3 , 



can be solved in the limit // 3> T, Ai, A2, when E, Ai and 
A2 can be considered homogeneous. In this limit, which 
is compatible with the self-consistent Born approxima- 
tion, the Green function (JjjJ is diagonal in frequency and 
momentum space and can be explicitly inverted: 

_ i0J n + £, p T 3 <7 3 + A1T1 + A 2 r 3 (Ti 



^ + ^ + A? + A| 

Here, we have defined the 'renormalized' expressions for 
the frequency and order parameters: 



1 - 



Ax 1 



A 2 1 + - 



r 




1 




2 ^ 


hi - 


FA 2 


+ A| 


r 




1 




2 ^ 




FA 2 


+ A| 


r 




1 




2 


M- 


FA? 


+ A 2 



A! 



= A 2 



(6) 



From the above equations of motion, one may deduce 
that 



Ax 



Ai ' 



(7) 



or, in other words, that in the weak disorder limit non 
magnetic impurities do not suppress s-wave supercon- 
ductivity (Anderson theorem^) while, introducing the 
parameters u = £D n /Aa and C = r/A 2 , 



A 2 



c 



v/l + u^l + A?/^) 



(8) 



Finally, the self-consistency equations for the SC and 
EI order parameters can be found minimizing the ac- 
tion J3J with respect to Ai j2 : 



7tAi ;2 \ - 



Ai, 



A? + A| 



(9) 



Here, Ai i2 = <?i, 2 ^f represent dimensionless coupling con- 
stants. Note that, as in conventional BCS theory, the in- 
tegral over momentum can be performed by making use 



of the identity / dp/(27r) 3 = / d£v(£,) ~ v F J d£. Em- 
ploying Eq. iJSJ, the self-consistency equations can then 
be rewritten in the form 



1 



E 



A? 



A 2 7T X - 



1 + V? 1 



A? 



-1/2 



-1/2 



(10) 



Combining together equations 10 with (|10fl . we are now 
able to discuss the finite and zero temperature mean-field 
phase diagram associated with the model JJJ. 



IV. PHASE DIAGRAM 

A. Temperature versus disorder phase diagram 

In the absence of disorder (i.e. T = and u — w„/A 2 ), 
one may note that, except for different coupling con- 
stants, the two self-consistency equations 111 01) are identi- 
cal. Therefore, since they cannot be satisfied simultane- 
ously, even though the SC and EI instabilities can occur 
simultaneously, in clean materials the corresponding or- 
derings are mutually exclusive. For Ai > A 2 the system 
becomes superconducting below 



7i(r = 0) = — 2O e 

7T 



-1/Ai 



(11) 



where £Iq is the frequency cutoff and — 1-78 is the 
Euler constant while, for A 2 > Ai, the transition into the 
EI state occurs at 



T 2 (r = 0) = — 2ft e~ 1/A2 



(12) 



Since charged non-magnetic impurities act as electron- 
hole pair breaking perturbations, while they do not af- 
fect the SC state, for Ai > A 2 the SC state dominates 
at any disorder strength T and the EI state never ap- 
pears. On the other hand, for A 2 > Ai, the EI phase is 
energetically more favorable at weak disorder, becomes 
suppressed for larger values of T, and eventually gives 
way to superconductivity. Non-magnetic impurities sup- 
press the EI state in exactly the same way as magnetic 
impurities suppress the SC statei 23 i 24 Therefore, one can 
infer that the dependence of the EI transition temper- 
ature T 2 on the disorder strength T is described by the 
standard Abrikosov-Gor'kov expression: 



In 



T 2 (T) \2 



27rT 2 (r)y V2 



(13) 



At some critical disorder strength r», the EI and SC tran- 
sition temperatures eventually coincide T 2 (r„) — Ti(T — 
0) = T* and, for r > T*, the system becomes supercon- 
ducting at the (r-independent) temperature T\ given by 
Eq. (TTJ. 
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FIG. 3: The temperature vs disorder strength phase diagram 
of the model Eq.Q for A^ 1 — AJ 1 = 0.5 . The dimensionless 
disorder strength is r/A 2 (T = 0, T = 0) and temperature is 
measured in units of T C 2{Y = 0). The EI and SC state are sep- 
arated by a thin domain, in which the two order parameters 
coexist. 



whereupon Eqs. (|10fl take the form 

1 1 ( 15 ) 

a 2 j^ v^+Af + r 

Note that the EI order parameter A 2 appears at tempera- 
tures lower than the "upper" T 2 (r) given by Ea. (|13fl . and 
disappears below the "lower" X^r), given by Eq. (TSJ. 
Figure 01 shows that Ti (at which the SC ordering sets 
in at T — 0) and r 2 (at which the EI ordering is de- 
stroyed at T = 0) are smaller than the disorder strength 
at the tetracritical point IV Therefore, in the interval 
T 2 < T < r„, the system passes through three con- 
secutive phase transitions as the temperature decreases: 
Firstly the system becomes an excitonic insulator, then 
it enters the mixed phase with the two coexisting order 
parameters, and finally the growth of the SC order pa- 
rameter with decreasing temperature suppresses the EI 
ordering, resulting in the transition into the pure SC state 
with A 2 = 0. 



One can therefore wonder in what way, at tempera- 
tures lower than T*, the transition between the EI and 
SC states takes place. In Fig. the temperature versus 
disorder phase diagram is shown for values of the cou- 
pling constants such that r\ = A7 1 — A^ 1 = 0.5. The pure 
EI and SC states are separated by a very thin region lo- 
cated in the T < T» and T < T* region of the phase 
diagram, where the two order parameters coexist. The 
three ordered phases (EI, SC, and EI+SC) and the high- 
temperature disordered phase merge at the tetracritical 
point (r»,T»). The boundaries between the coexistence 
region and the two pure phases are critical lines of second 
order transitions although, due to the very small width 
of the coexistence region, the evolution of one pure phase 
into another is close to being of first-order. This will be 
further discussed in Sec. 

At the first critical line Ti(r) separating the pure EI 
phase from the mixed EI+SC phase, one has Ai — > 
and A 2 ^ 0. In the Ai — > limit, the self-consistency 
equations (|10[) simplify to the expressions 



1 „ \ -> u 
— = 2ttTi V == 

^ = ^£^===. 



(14) 



The second critical line r 2 (r) separates the mixed state 
from the pure SC state and is obtained by instead taking 
the limit A 2 — > for nonzero Ai. In this case, one can 
approximate 




B. The zero temperature phase diagram 

The zero temperature phase diagram is shown in Fig.0] 
The EI state exists only for positive ry (or equivalently 
e-'' = Ai(T = 0)/A 2 (T = 0, T = 0) < 1). The coexis- 
tence region (shaded) is confined between the two criti- 
cal lines Ti(rj) < T 2 (r]). The system is superconducting 
for r > T\, while the excitonic condensate appears for 
r < r 2 ). For small 77, i.e. close to the quantum critical 
point separating the EI and SC states in the absence of 
disorder, the critical disorder values can be obtained as 



72(77) 



-77 



-77 



3^+8 2 
3ir 3 'I ' 



(16) 



where 7 = r/A 2 (0,0). Therefore, the width of the coex- 
istence region is approximately given by 7 2 — 71 ~ 0.08ri 2 . 

For large values of 77, i.e. when the SC coupling Ai 
is much smaller than the EI coupling A 2 , the coexis- 
tence region essentially coincides with the disorder inter- 
val e _7r//4 < 7 < 1/2 in which the EI state is gaplessi 23 i 24 
Therefore, for 77 > 1 + 37r/4, the superconductivity ap- 
pears at the same disorder value, at which the EI becomes 
gapless, 71 = e _7r / 4 ~ 0.46, while T 2 asymptotically ap- 
proaches the disorder strength at which the EI state is 
destroyed in the absence of superconductivity: 



7 2 ~ 1/2 - rye' 



-2v 2 



i] > 1 



(17) 



We note that, for rj < l+37r/4, the gap in the spectrum of 
quasi-particle excitations at zero temperature is nonzero 
for all r, while for 77 > 1 + Sir/ 4 it vanishes at a single 
point r = T±. 
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FIG. 4: The zero temperature phase diagram. Close to the 
quantum critical point r\ — V — 0, the coexistence region is 
extremely thin. For i/ < 1 it practically coincides with the 
disorder interval, in which the EI is gapless. 



This re-entrant behavior and the form of the phase di- 
agram are similar to what was found for the spin-Peierls 
compound CuGe03 , which upon doping shows an anti- 
ferromagnetic ordering coexisting with spin-Peierls phase 
in some interval of doping concentrations^^ 



V. EI-SC SYMMETRY 

To understand why the coexistence region is so narrow, 
it is instructive to plot the EI and SC order parameters 
as functions of F in the coexistence region (see Fig. 03). In 
the interval ri < T < r 2 the excitonic (superconducting) 
order parameter decreases (increases) fast with increasing 
r, while A = tj A| + A 2 , stays approximately constant. 

This behavior results from the symmetry between the 
6- and a-electrons at the quantum critical point Ai = A2 
in the absence of disorder (r = 0) : 



(18) 



the replica trick: 



47T x - 



A? 



2vrr 



A 2 
^2 



A? 



A 2 
^2 



A 2 
^2 



A 



sign(wn)- 



A, 



"An 



(19) 



in the absence of disorder , th e free 
VAf 



Let us notice that 

energy only depends on 'the total gap' A 



A 2 

^2- 



This follows from the fact that the generator of the EI-SC 
rotations T2U1 commutes with the Hamiltonian Eq. © 
for U(x) = 0. Moreover, for g\ — gi, the last term in 
Eq. l(T9l is equal to 2A 2 /gi (the second term in Eq. JT3J 
vanishes for T = 0). Thus, at the quantum critical point 
the free energy has a 'Mexican hat' profile as a function 
of the order parameters (Ai,A2), symmetric under the 
rotations transforming the excitonic insulator into the su- 
perconductor. This symmetry between electron-electron 
and electron-hole pairing is analogous to the symmetry 
unifying the d-wave superconductivity and antiferromag- 
netism discussed in the context of high-T c and heavy 
fermion materials. 29 ' 30,31 

Away from the quantum critical point, and for nonzero 
disorder, the electron-hole symmetry is broken. Solving 
Eqs. © for tint Ai, and A 2 , perturbatively in the disor- 
der strength F, and replacing at T = the summations 
over the Matsubara frequency uj n in Eq. I|19fl by inte- 
grals, one can obtain an expansion of the average energy 
density in powers of T: 



Ml 

2v F 



A 2 A 



Ai 



A, 



Mo 



^A 2 1A 2 (3A 2 + A 2 ) 



h -r — 

2 A 3 t±- 
n A 2 A 2 (A 2 -4A 2 ) . : 



r- 



16 



A 7 



r 3 + o(r 4 ) . (20) 



Denoting the dimensionless disorder strength by S 
r/ A -C 1 and defining the angle cf> 

Ai = Acos0 A2 = Asint/) , 

we can recast Eq. (I2UII in the form 



CO 

2v F 



A 2 



In ^ - - + - cos 20 + — (1 - cos 20) 
A 2 2 4 v ; 



This transformation results is the rotation in the space 
of the two order parameters over the angle <fi E [0, 2ir], 

Ai 1— > Ai cos <f) — A2 sin 4> 
A2 1— > A2 cos + Ai sin . 

At the mean-field level, the anomalous part of the av- 
erage free energy per unity of volume, (/) = (F)/V, can 
be easily evaluated starting from and making use of 



— — (3 — 2 cos 20 — cos ' 



tt6 3 

+ — — (6 + 5 cos 20 — 6 cos' 



5cos60) + O((5 4 



(21) 



where A is the geometric mean of the EI and SC or- 
der parameters, A = N /A 1 A 2 (0,0) = 2fi e"( A ^ 1 + A ^ 1 )/ 2 . 
The symmetry-breaking terms in Eq. I|21|l (that depend 
on the angle 0) are proportional to powers of r\ and S. 
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FIG. 5: The plots shows the dependence of the superconduct- 
ing order parameter Ai (dashed line) and the excitonic order 
parameter A2 (thin solid line) on the disorder strength V in 
the region of coexistence of the two phases for 77 = 0.1 and 
T = 0. For Ti < T < T 2 , Ai and A 2 vary very fast, while 
A = \J Ai 2 + A2 2 (thick line) stays approximately constant. 
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FIG. 6: The plot shows the energy minimum versus the angle 
(j> at r\ — 0.5 and for different values of V: T = 0.234 < Ti 
(solid line), T ~ F fo = 0.237 (dashed line) and F = 0.240 > 
T2 (dot-dashed line). 

Thus, for 77, S <C 1, these terms are small and the energy 
has the slightly deformed 'Mexican hat' shape with an al- 
most flat valley connecting the points A m ;„(0), at which 
(/) has a minimum for a given <fi. 

This is illustrated in Fig. where we plot the <f>- 
dependence of the minimal energy density for, respec- 
tively, T = 0.234 < Ti, T = 0.237 (in the coexistence 
region) and V = 0.240 > T 2 (in both cases 77 — 0.5). 
For r < Ti, the energy minimum is at 4> — 7r/2 (the EI 
state), while at T > r 2 the energy minimum is for = 



(the SC state). Though the ^-dependence of the minimal 
energy is, in general, rather complicated, the scale of the 
energy variations in all three cases is very small, i.e. the 
valley is practically flat. This is the reason for the narrow 
width of the disorder interval ri < T < r 2 , in which the 
two phases coexist - a very small variation of the disor- 
der strength 7 is sufficient to shift the position the energy 
minimum from </> = 7r/2to0 = O along the energy valley, 
in which A = yj Af + A| remains practically unchanged. 

FigureEJalso illustrates the absence of first-order tran- 
sitions in the modelQ). Note that T = 0.237 is close 
to Ffo, at which the energies of the SC and EI states 
become equal: (/(0)) = (/(7r/2)). The first-order tran- 
sition between the two pure states, however, does not 
occur, since the energy has the global minimum at some 
angle <fi, such that < (f> < tt/2, corresponding to a mixed 
ground state. Furthermore, when T — 0.237, the energy 
has a local maximum at = tt/2, enforcing the EI state 
to be metastable. 



VI. DISCUSSION AND CONCLUSIONS 

We discussed effects of disorder in systems with com- 
peting instabilities, such as CDW superconductors. We 
considered a simple model, which describes a metal with 
two perfectly nested electron-like and hole-like parts of 
the Fermi surface. In this model the interplay between 
the electron-electron and electron-hole pairings is very 
strong, as these they compete over the whole Fermi sur- 
face. 

We showed that disorder can be used to tune the bal- 
ance between the two competing phases and to stabilize 
the state, in which they coexist. The charged nonmag- 
netic impurities induce superconductivity by suppressing 
the CDW state. Such a disorder-induced superconduc- 
tivity is observed in the irradiated two-dimensional CDW 
material 2ff-TaS 2 »i2i In other transition metal dichalco- 
genides, e.g. 2if-NbSe 2 and 2iJ-TaSe 2 , which are CDW 
superconductors already in absence of disorder, a small 
amount of irradiation-induced defects results in an en- 
hancement of T c m^ Similar behavior is observed in the 
quasi-one-dimensional CDW material Nb 1 _. E Ta. l; Se3. In 
the pure NbSe3, T c is smaller than 50 mK at ambient 
pressured The substitution of Nb for Ta suppresses the 
resistivity anomalies due to the CDW transitions, while 
T c grows up to ~ 2K at x — 0.05£ The effect of impuri- 
ties in these materials is similar to that of pressure and 
hydrogen intercalation. 5,16,17 

In agreement with these experimental findings, the 
phase diagrams of our model Figs. and show a strong 
sensitivity of the ground state to disorder and the coexis- 
tence of the SC and EI states in the presence of disorder. 
This behavior can be easily understood and described an- 
alytically, using the Landau expansion of the free energy 
in powers of the EI and SC order parameters near the 
quantum critical point [see Ea. (|20|) ]. which we derived 
from the microscopic model. Disorder distorts the shape 
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of the energy potential and continuously shifts the posi- 
tion of the minimum from the point corresponding to the 
excitonic insulator to the point corresponding to the su- 
perconducting state, which gives rise to the coexistence 
of the two states. 

The microscopic origin of this coexistence is the break- 
up of a part of the electron-hole pairs by disorder and the 
subsequent recombination of the released fermions into 
electron-electron and hole-hole pairs. In other words, 
disorder transforms the CDW gap in the single-electron 
density of states into a pseudogap, filled with states de- 
scribing the broken electron-hole pairs. The SC phase 
develops inside this pseudogap, which resembles the be- 
havior observed in high-T c cuprates. 32 

In addition to the disorder-induced superconductivity, 
resulting from the suppression of the EI state, the phase 
diagram of our model (see Fig. shows an interesting 
'inverse' effect, namely the suppression of the EI state 
due to the growth of the SC order parameter with de- 
creasing temperature. Though this re-entrance transi- 
tion is just another consequence of the competition be- 
tween the two types of ordering, we did not find any 
reports of such a behavior in CDW superconductors in 
the literature. This, however, it has been observed in the 
quasi-one-dimensional spin-Peierls compound CuGe03, 
where impurities induce the long range Neel ordering*^ 
In this material, the interplay between the dimerized and 
antifcrromagnetic states allows for a similar theoretical 
description^ In most CDW superconductors the CDW 
transition occurs at a much higher temperature than the 
SC transition, so that the influence of the Cooper pair- 
ing on the CDW modulations is difficult to observe. Fur- 
thermore, the CDW gap only opens on a nested part of 
the Fermi surface. In quasi-one-dimensional NbSe3 the 
fraction of the Fermi surface affected by the CDW tran- 
sition was estimated to be ~ 0.6 at ambient pressure.^ 
In the two-dimensional 2i7-NbSe2 this fraction is appar- 
ently very small, since the gap opening actually increases 
the conductivity of this material^ and the part of the 
Fermi surface, where the gap opens, was not found in 
ARPES experimentsji 3 ^^ even though the gap value 
(34meV) is knowniM For partially gapped Fermi surfaces 
the competition between the CDW and SC states is less 
strong, so that in 2_H-NbSe2 they coexist even in absence 
of disorder. 

While the enhancement of the SC transition tempera- 
ture upon the suppression of the CDW state is well doc- 
umented in many materials, the experimental situation 
with influence of the superconductivity on the CDW state 
is less clear. On the one hand, Raman experiments on 
2ii-NbSe2 show the suppression of the intensity the col- 
lective SC mode by magnetic field with the concomitant 
enhancement of the intensity of the CDW modesiSiSi On 
the other hand, no effect of the superconducting ordering 
below 7K and of the suppression of the SC state by mag- 
netic field on the CDW modulation was observed in x-ray 
experiments. 1 ^ The understanding of the behavior of 2H- 
NbSe2 is complicated by the multi-sheet structure of the 



Fermi surface and the momentum- and sheet-dependence 
of both order parameters^ The interplay between the 
CDW and SC states in this and other materials requires 
further experimental and theoretical studies. 

Crucially, one may note that the phase diagram of the 
interacting system was inferred from the self-consistent 
Hartree-Fock approximation which captures only the 
mean-field characteristics. In view of the filamentary 
structure of the coexistence region, the system can be 
susceptible to mesoscopic or sample to sample fluctua- 
tions due to the quenched impurity potential. Such ef- 
fects are recorded in fluctuations of the field £(x) around 
its saddle-point or mean-field value (as opposed to the 
leading terms gathered in the low-order T expansion con- 
sidered here) . In the gapless regime, such effects can give 
rise to long-ranged diffusion mode contributions to the 
generalized pair susceptibility (cf., e.g., Ref. How- 
ever, in the present case, the disorder potential imposes 
a symmetry breaking perturbation on the EI phase. As 
such, we can expect mesoscopic fluctuations due to dis- 
order to impose only a short-ranged (i.e. local on the 
scale of the coherence length of the EI order parameter) 
perturbation on the pair susceptibility. In the vicinity of 
the coexistence region, where the potential for the angle 
4> is shallow, the effect of these mesoscopic fluctuations 
may be significant. 

To understand the effect of random fluctuations in the 
coexistence region, we consider the Ginzburg-Landau ex- 
pansion for the ground state energy close to the quantum 
critical point 77 = 0, T = 0, at which the EI and SC states 
are degenerate. In the vicinity of this point the phase <fi 
of the 'total' order parameter Ai + 1/S.2 = Ae 1 ^ is a soft 
mode, so that weak disorder mainly induces spatial fluc- 
tuations of the phase, while the magnitude of the order 
parameter A approximately stays constant. As in the 
derivation of Eq. 12I|) , we expand the energy in powers of 
A and disorder strength, assuming that r <C A, which is 
justified in the coexistence region, where T ~ ^A (see 
Ea. (|I6|) . Assuming that the phase varies slowly on the 
length scale of the correlation length £ = ^ , where vf 
is the Fermi velocity, we obtain 



F : 



f(Vor • ,/:x;ms2o) (/) 



(22) 



where the first term describes the 'elastic energy' of an 
inhomogeneous state, the disorder-averaged free energy 
(/) is given by Ea. (|2*T)l . and u(x) is the fluctuating part 
of disorder coupled to the phase of the order parameter. 
Neglecting correlations on a scale smaller that the corre- 
lation length, u(x) can be approximately considered as a 
random 5-correlated Gaussian variable with zero average, 
(m(x)) = 0, and variance 



(u(x)u(y)) =AS(x-y) 



A = 



tt 4 T 2 A 2 
2kp 



(23) 



(we omit the lengthy calculations that lead to this result). 
The coupling to disorder also occurs in higher orders of 
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the expansion, but those terms are relatively small and 
can be neglected. 

Following the Imry and Ma argument^ we consider a 
large phase fluctuation, e.g., a droplet of the SC phase 
of the spatial extent L inside the EI matrix. Comparing 
the ty pical energy gain due to the coupling to disorder 
~ vp v AL 3 with the loss in the elastic energy ~ i/pv F L, 
we find that the fluctuation is energetically favorable for 



L > 



Vp_ 

A 



k 3 



Jp^p 



27T 2 77 2 A 4 



(24) 



where we took into account that, in the coexistence re- 
gion, r ~ ^A. _ 

The crucial difference of our model from that consi- 
dered in Ref.[38| is the absence of an exact continuous 
symmetry. Even in the coexistence region, the minimal- 
energy valley connecting the SC and EI points ((f) — 
and 4> — § ) is not perfectly flat. The typical amplitude 
of the variations of the energy density is ~ rfvpA 2 (see 
Eq. J2U), resulting in the energy loss ~ r] 2 vpA 2 L 3 pro- 
portional the volume of the fluctuation, which suppresses 
large droplets. Comparing it with the energy gain, we 
find 



L A < 



A 



?7 4 A 4 



Equations l|24[l and (|25[1 hold simultaneously for 

Vpkp < T]A , 



(25) 



(26) 



which cannot be satisfied in the weak coupling limit. One 
may wonder why the condition l|26|l does not hold even 
for rj = r = 0, where the model has a continuous symme- 
try. The reason is that in our model the role of disorder is 
two-fold. On the one hand, it couples to the order param- 
eter, as in the 'random field' model discussed in Ref.pBSl] 
and tends to destroy the ordering. On the other hand, 
it affects the energy difference between the EI and SC 
states and, therefore, suppresses the phase fluctuations, 
by destroying the symmetry of the energy potential. The 
second effect, which is linear in T, is stronger than the 
first. 

Thus, the inhomogeneity of the order parameter, re- 
sulting from typical disorder fluctuations is small. The 
phase fluctuations can also be induced by large disor- 
der fluctuations ('Lifshitz tails'), but their contribution 
to the free energy is exponentially smalli^ 9 . This justifies 
our mean field treatment of disorder. 

This conclusion may not hold, however, for strongly 
coupled CDW superconductors or for other types of dis- 
order. Qualitatively, we expect that inhomogeneous ex- 
citonic and superconducting order parameters may result 
in a broadening of the coexistence region. The local sup- 
pression of the excitonic pairing near charged impurities 
can give rise to the local enhancement of the supercon- 
ducting order. The state with such a nanoscale phase 
separation, in which two competing orders alternate in 



antiphase without a loss of the macroscopic coherence, 
can be more energetically favorable than the uniform 
state and, therefore, can be stabilized in a wider interval 
of parameters. Such a state was observed in /iSR experi- 
ments on doped CuGe03, which shows both spin-Peierls 
and antiferromagnetic ordering. 40 

In conclusion, we studied effects of disorder on systems 
with competing superconducting and chargc-density- 
wave instabilities. We showed that even in the extreme 
situation, when the competition takes place over the 
whole Fermi surface and the superconducting and charge- 
density-wave phases are mutually exclusive, disorder can 
give rise to their coexistence in a spatially homogeneous 
state. Furthermore, disorder itself can be used as a pa- 
rameter, with which one can tune the balance between 
competing phases. Although our model is too simple 
to describe the physics behind the coexistence of super- 
conductivity and CDW(SDW) states in, e.g., high-T c or 
heavy fermion materials, we believe that the ability of 
disorder to bring together incompatible phases may be 
important for understanding phase diagrams of these sys- 
tems. 



APPENDIX A: COEXISTING INSTABILITIES 
AND DERIVATION OF THE EFFECTIVE 
MODEL 



In this appendix we obtain a condition under which 
both the SC and EI instabilities can occur simulta- 
neously. Here, we consider more realistic interac- 
tions between electrons than those described by the 
model l[T|l. namely, the phonon- mediated interaction and 
the Coulomb repulsion. The Coulomb repulsion counter- 
acts the phonon-mediated attraction between electrons 
and suppresses the SC instability. The same holds for the 
instability towards the formation of the excitonic conden- 
sate with the difference that the two types of interaction 
now change roles: the Coulomb force favors the electron- 
hole pairing, while the one-phonon exchange results in a 
repulsion between electrons and holes. We will show that 
the SC and EI instabilities can coexist due to different 
frequency dependence of the two types of interactions. 

For retarded phonon-mediated interactions, the order 
parameters Ai and A2 are frequency-dependent, which 
complicates the solution of the self-consistency equations. 
We show, however, that in the weak coupling and weak 
disorder limit, the equations for the order parameters at 
zero frequency coincide with Eq. ©, which justifies the 
model introduced in Sec. [H] Moreover, we will give the 
explicit expressions for the coupling constants g± and (72 
appearing in Eq. JQ. 

We describe effective electron-electron interactions by 
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a non-local action 

S cS =^jdxj drdT'p( X ,T)D(r~T')p( X ,T') 

+ 9hJ dxdrp 2 (x,r) , (Al) 

where p = X^OLVw + ip ba ipba) is the total electron 
density. The first term is the phonon-mediated effective 
attraction between electrons and D{t — t') is the phonon 
Green function. For a single dispersionless optical 
phonon with the frequency f^o and the propagator D u!n — 
-ng/(w2 + n 2 ), we have D(r - r') = -ft e~ 0o|r ~ T,| /2 
for T <C SIq. The second term in Eq. ijAlfl is the in- 
stantaneous Coulomb interaction. We neglect the mo- 
mentum dependence of the screened electron-phonon and 
Coulomb couplings, which makes the electron-electron in- 
teractions local in space. 

The couplings for the a and 6-electrons in Eq. IjAlfl give 
rise to a large freedom in the choice of order parameters, 
which in reality may not be present, e.g., due to the inter- 
band scattering, which separately does not conserve the 



numbers of the a— and b— electrons. In what follows we 
restrict ourselves to the anomalous averages considered 
in Sec.[n] which, for retarded interactions Eq. IjAljl . are 
time-dependent (j = a,b and a =|, J.): 

Ai(r -r') = - [gi ph D( T - r') + g 2 c ] (^ T (r)V a (r')> 
A 2 (r - r') = - [gi ph D(r - r) + g 2 c ] (^(r)^(r')) . 

In the frequency representation the self-consistency equa- 
tions read 

Ai(w„) = - — (5o- P h^„-< - 9c) tr (tiG Pi< ) 
A 2 (w n ) = — Y (3c-ph-CU-< - 9c) tr (r 3 CTiG Pi< ) , 

where the electron Green function is given by Eq. (J^J . 

To simplify the algebra, we consider here only the zero 
temperature case. The integration over the electron ex- 
citation energy £ gives 



1 r +Ee 
AiM = - J duo' 



1 f +E " 
Ai(w) = - / dJ 

1 J-E a 



Kl 



n 2 



(uj - uj') 2 + n 2 



K 2 



n 2 



{uj - uj 1 ) 2 + n 2 



Cu' 2 + A 2 (uj>) + A 2 (uj') 
uj' 2 + A 2 (uj') + A 2 {uj') 



(A2) 



where we have introduced the dimcnsionless coupling 
constants, K\ = ^F<7o-ph an< ^ K 2 = v f9ci an( i where 
E c is the frequency cutoff required for the instantaneous 
Coulomb interaction. Moreover the variables Q,Ax, and 
A2 are defined in @. 

Although Eqs. (|A2|) look at a first sight complicated, 
one can see that, in the limit of weak coupling and weak 
disorder, Ai l2 ,r <C fio < E c , their solution can be found 
by making use of the fact that the order parameters 
Ai(w) and A 2 (w) strongly vary at frequencies uj ~ flo, 
while uj, Ai, and A2 are nontrivial functions of uj only at 
much lower frequencies uj ~ T, where A\(lu) and A 2 (uj) 
can be replaced by their zero frequency values. There- 
fore, we can solve Eqs. <|A2ll in two steps: first we find 
the frequency dependence of the order parameters Ai(cj) 
and A 2 (uj) for arbitrary values of Ai(0) and A2(0), and 
then we solve the self-consistency equations for Ai(0) 
and A 2 (0). 

It is convenient to use the dimensionless variables x — 



uj/Qq and 2/1.2 (x) = Ai.2(w)/f2o (and similarly x and 
2/1,2(2;)), in terms of which the first of the equations (|A2|) 
reads 



j/i (x) = / dx' 



yi{x') 
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2 



^/x' 2 +yl (x') +y%(x') 
1 1 



(x + x') 2 + 1 (x-x') 2 + l 



K 2 



where A = E c /Qq. We then introduce an intermedi- 
ate scale X, such that 2/1,2 *C X <C 1. In the interval 
< x' < X, we can neglect the x'-dependence of the 
kernel of this integral equation and the functions 2/1,2(2;') 
(however, j/1,2 and x' still do depend on x'). In the second 

interval X < x' < A, we substitute 2/1/ y x' 2 +2/1+2/2 by 
yi(x')/x' and perform the integration by parts. In this 
way we obtain 
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Vi( x ) = -2/1 (A) «2 In A + 



«i 



K 2 



da;' 



yf (a;') + 2/2V) 



2/1 (0) In J*f 



da/ In : 



dx' 



Kl 

2 



(.x + x') 2 + 1 {x - x') 2 + 1 



- K 2 



2/1 (»') , (A3) 



where the limits of the second integration were extended 
to and oo, as there is convergence both at small and 
large frequencies. 

Since, at X ^> yi i2 , Eq. i|A3fl is independent of X, we 
can chose X — 1 (and still substitute ?/i j2 (x') by 2/1,2(0) 
in the first integral) . The value of yi at the cutoff is then 
given by 



2/i (A) 



dX ; 



Si(s') 



yf(x') + 2/2V) 



dx In x 



d_2/i 
dx 



where = k 2 /(1 + /c 2 In A). For arbitrary x we have 

Kl 



yi{x) 



1 



«2 %i(Q), 2/2(0)) 



dx' In x' - 



dx' 
1 



Kl 

2 



1 



(X-X') 2 + 1 



(x + x') 2 + 1 

2/1 O^O 



(A4) 



where the notation 



m(x') 



7(2/1(0), 2/2(0))= / dx' — , , ^_ 

do v/x' 2 + y 2 (x') +^(x') 

is used to stress the fact that y\ and y 2 are assumed to 
be frequency independent. 

In the weak coupling limit the first term in the right- 
hand side of Eq. (|A4|) . proportional to the 'large loga- 
rithm' In \/y\ + y 2 , is much larger than the second term, 
so this integral equation can be solved by iterations, 
which generate a perturbative expansion for yi(x). To 
the lowest order, the frequency dependence of the order 
parameter coincides with that of the kernel* 2 ^ 



2/1 (x) = 



J(tfi(0) >ttl (0)) 



(A5) 



Then the self-consistency equation for Ai(0) coincides 
with Eq. @ at T = 0: 



Ai(0) = Ai / duo 
Jo 



/J) 2 + A 2 (^) + A 2 (w) 
and the effective coupling constant is given by 

K 2 



Al = Kl 



Kl 



1 + k 2 ln_E c /51o 



(A6) 



(A7) 



The negative term in the coupling constant describes the 
reduction of the attraction between electrons due to the 
Coulomb repulsion, but this reduction is itself reduced 
by the presence of the large logarithm in the denom- 
inator due to the difference in the time scales of the 
retarded phonon-mediated attraction and the Coulomb 
rcpulsion^ii The first-order correction to yi(x), found 
by substituting Eq. (|A5|) into the integral in Eq. I|A4|) 
(as well as all higher-order corrections), leaves the form 
of the self-consistency equation (|A6|) unchanged, but re- 
sults in a small modification of the expression for the 
effective coupling constant through Ki and /c 2 : 



Al = Kl 



The frequency dependence of the excitonic insulator 
order parameter A 2 (u>) and the self-consistency equation 
for A 2 (0) can be obtained from Eas. lA6IA7ll by the sub- 
stitution Ai 1— > A 2 , ki 1 — ► — Ki, and k 2 1— > — k 2 : 



mo) = a 2 r du 

J " ^ 2 + a 2 h + a 2 h 



(A8) 



where the effective coupling constant A 2 is to the lowest 
order given by 



A? 



K-2 



1 — K2 In E c /Qq 



Kl . 



(A9) 



In Fig. [21 we show the typical frequency dependence of 
the SC and EI order parameters, calculated for Ki = 0.25 
and Ki = 0.1 and T = 0. (Since in the absence of disorder 
the SC and EI states cannot coexist, we calculated Ai (oj) 
assuming A 2 = and vice versa.) The SC order param- 
eter Ai(ix>) is positive at small frequencies and changes 
sign at u> ~ fio 5 while the EI order parameter A 2 (w) has 
a dip for \uj\ < FIq. The 'separation' of the two order 
parameters in frequency is crucial for the coexistence of 
instabilities. 

The necessary condition for superconductivity to ap- 
pear is Ai > 0, while the instability towards the excitonic 
condensate occurs for A 2 > 0. These two conditions, 



Ai > 



A, 



1 + A 2 In A 
hold simultaneously for 



Ai < 



A, 



1 - A 2 In A 



r'/l 



1 



1 



<ln n 



(A10) 
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A weak disorder has little effect on the frequency de- 
pendence of Ai and A 2 . However, its presence is crucial 
for the stabilization of the mixed state, in which the two 
order parameters coexist. 
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